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A condensed introduction to the basic concepts of causal perturbation theory is given. Causal per- 
turbation theory is a mathematically rigorous approach to renormalization theory, which makes it 
possible to put the theoretical setup of perturbative quantum field theory on a sound mathematical 
basis by avoiding infinities from the outset. It goes back to a seminal work by Henri Epstein and 
Vladimir Jurko Glaser published in 1973, where a specific causality condition was imposed at 
every order of perturbation theory in the case of scalar quantum field theory such that divergent 
integrals could be avoided in actual calculations of loop diagrams. In the meantime, the causal 
approach has been applied also to a wide range of gauge theories. 
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1. Introduction 

Before we come to the actual problem of ultraviolet (UV) divergences in perturbative quan- 
tum field theory (pQFT), we discuss here a naive example of 'UV divergence' by considering 
Heaviside-0- and Dirac-5-distributions in 1-dim. 'configuration space'. The product of these two 
distributions 0(jc)5(jc) is obviously ill-defined, however, considering the Fourier transforms 

&{8}{k) = S(k) = J dx8{x) e - ikx = 1, (1.1) 



&(k) = lim / dx&(x)e- ikx - ex = lim 

e\0j e\0 k-lE 



k-iO 
o 



(1.2) 



one can nevertheless calculate in a formal manner the Fourier transform of the ill-defined product 
mentioned above 

3?{eS}{k) = J dxe- ikx e{x)S{x) = J dxe- ikx J ^@{k')e +ik ' x J ^8{k")e +ik " x . (1.3) 
Since / dxe'^ k ' +k "^ x = 2n8{k' + k" — we obtain the divergent convolution integral 

*{«»<*) -^fde^k-n) = -i-J^- v (i.4) 

The obvious problem in x-space leads to a divergent integral in k-space. This situation arises in a 
completely analogous manner in pQFT, as will be discussed in the following. 



2. The origin of UV divergences in perturbative quantum field theory 

In pQFT, the role of the Heaviside ©-distribution is taken over by the time-ordering operator. 
The well-known textbook expression for the perturbative scattering matrix given by 

-[-co ~\~°° 

S=f, I dh... I dt n T[H int (h) . . .H int {t n )\ 

— CO —CO 

= £ ^T- I d ^ ■ ■ ■ I d 4 x n T[J% m ( Xl ) . ..J% nt {x n )l (2.1) 

where the interaction Hamiltonian Hi nt {t) is given by the interaction Hamiltonian density J#i nt (x) 
via Hi nt (t) = j d 3 xJ^ nt (x), is problematic in the UV regime (and in the infrared regime, when 
massless fields are involved). A time-ordered expression a la 

T[J% m (x 1 )...J4? int (x n )}= £ 0(xi\ 1 -xi ) l2 )...0(x° (n _ i) -xl , I J^(xn 1 )...^(xn„) (2.2) 

Perm. Yl 

is formal (i.e., ill-defined), since the operator-valued distribution products of the are simply 
too singular to be multiplied by ©-distributions. This observation provides a formal explanation 
for the existence of UV divergent expressions in pQFT, which are usually related in a qualitative 
manner to contributions of virtual particles with 'very high energy', or, equivalently, to physical 
phenomena at very short distances. 
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A first step towards the solution of the apparent mathematical problem was taken by N. N. Bo- 
goliubov and D. V. Shirkov [1], which introduced a clear definition of the causality condition in 
pQFT and the concept of adiabatic switching (see below). However, in their attempt the UV diver- 
gences did persist. A rigorous mathematical analysis within the framework of distribution theory 
was finally presented by H. Epstein and V. Glaser [2], who derived an inductive construction of 
the perturbation series based on Poincare invariance and causality (unitarity plays no immediate 
role). In their approach, UV divergences are avoided from the start, and the Feynman rules only 
hold on tree-level. The calculation of loop diagrams turns out to be rather technical and involves 
finite (subtracted) dispersion integrals instead of divergent Feynman integrals. A new strategy to 
treat the infrared problem by adiabatic switching of the interaction is also introduced in the causal 
approach. 

In the following, an introduction to the basic concepts of the causal Epstein-Glaser method is 
given. It is important to note that the presentation below suffers from a certain lack of mathematical 
rigor as a natural consequence of the limited space available in these proceedings. For a concise 
introduction to the causal method we refer to the textbook of G. Scharf [3]. 

3. Mathematical preliminaries 

3.1 Free quantum fields and operator valued distributions 

It is a crucial observation that free field operators are operator- valued distributions. E.g., for a 
scalar (neutral) field of a particle with mass m one has 



where the annihilation- and creation operators a and a + fulfill the usual commutation relations. 
This expression must be smeared out by rapidly decreasing test functions g(x) in the Schwartz 
space Y(R 4 ) (5? (R 3 ) would also be sufficient in the case of a free field) in order to get an operator 
in Fock space, formally written in integral form as (p(g) = J d 4 x(p(x)g(x). Note that (p(x)\0) is not a 
Fock state. The same arguing applies to the interaction Hamiltonian densities used in perturbation 
theory constructed from normally ordered products of free fields. E.g., in QED the perturbative 
interaction Hamiltonian density is expressed by the help of the free spinor field ^(x) and the photon 
field Afi{x) by the well-defined operator- valued distribution 3tf nt = — e : x ¥(x)j fix i'(x) : A^(x). In 
<p 3 -theory, which will serve as a model theory in the following, one has Jtfj nt = |y : (p(x) 3 :. 

3.2 The perturbative S-matrix: basic properties 

Based on the observations made above, it is therefore most natural to replace the problematic 
expression eq. (2.1) by 



where the switching function g(x) can be considered as a local variation of the coupling constant 
and we set T\(x) = —iJifl nt (x). In the causal approach, the T n (x\, . . . ,x n ) ~ T[T\(x\) . . .T\(x n )\ 




(3.1) 




n=l 



(3.2) 
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denote a well-defined (divergence free) time-ordered product, which is symmetric in the sense that 

T n \- . . ,xi, . . . ,xj, . . .) — 7Jj(. . . ,xy, . . . ,Xj, . . .) Vi, j (3.3) 

by construction. Since the T„ are constructed such that they are free of any UV divergences, every 
order of the perturbative S— matrix S n (g) = X J d A x\ . . .d 4 x„ T n (x\,... ,x n )g{x\) . . -g{x n ) is well- 
defined even when massless fields are present. Infrared divergences are absent as long as the 
interaction is switched by g G ^(B 4 ). Infrared problems arise in the so-called adiabatic limit 
g(x) — > 1, since 1 ^(IR 4 ). Performing the adiabatic limit is a delicate task as far as existence 
and uniqueness of the limit are concerned. The limit has to be taken such that observable quantities 
(like cross sections) remain finite. A typical strategy is to rescale g(x) according to 

limg(ex) "— g{x) =g(0) = const. (3.4) 

e\0 

No further regularizations, like the introduction of a finite photon in mass in theories like QED, are 
then necessary in this approach. 

One should note that no statements about the convergence of the full series eq. (3.2) can 
be made in general. We further mention that perturbative expansion of the inverse perturbative 
S-matrix is given by 

S(g)- 1 = 1 + £ - / d A x x . . . d\ f n (xi , . . . ,x n )g( Xl ) . ..g(x n ) = (t + T)- l =l + £ (-T) r , 

n=\ n - •' r=l 

(3.5) 

n 

where T n (X) = £ (— l) r L^«i C^i) • ■ ■ T nr (X r ), and X = {xi ,x„} is a disordered set and £ de- 

r=l P r P r 

notes all partitions of X into r disjoint subsets 

X=X l U...UX r , Xj^Q), X i f]X j = Q), \Xj\=iij. (3.6) 

4. The method of Epstein and Glaser 

4.1 Inductive construction of the perturbative S-matrix 

Causality is the pivotal point in the UV divergence-free approach of Epstein and Glaser. 
Causality is expressed by the condition that if the switching function g(x) = g\(x) + gz{x) can 
be decomposed such that the supports of g\{x) and g2{x) (denoted by supp(gi) and supp(g2)) are 
space-like separated, i.e. if there exists a reference frame such that from xGsupp(gi) follows 
x° < and fromy Gsupp^) follows y° > (see Fig. 1), then one has S(gi +#2) = S(g2)S(gi) for 
all switching functions gi,g2 which fulfill the condition above (denoted by supp(gi) < supp(g2), 
where the symbol '<' should be read as 'earlier'). 

This implies T n (xi, . . . ,x n ) = T m (xi,. . . ,x m )r„_ m (x m +i, . . . ,x„) if {x\, . . . ,x m } > {x }, a 

condition which is, of course, intuitively clear. 

Note that the support of a (test) function supp(g) is the closed set obtained by taking the 
complement of the largest open set on which g vanishes. Thus we adopt an analogous definition of 
the support for distributions, which is the complement of the largest open set on which a distribution 
vanishes. A distribution vanishes on an open set if it vanishes for all test functions whose supports 
are in the open set. 
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Fig. 1: An example for supp(g 2 ) > supp(gi). 
4.2 Explicit construction of 82(g) 

Causality and translation invariance require that the commutator D 2 (z = x\ —X2) 

£> 2 (x 1 -x 2 ) = (-/) 2 [^ f (x 1 ),^(x 2 )] = [r 1 (x 1 ),r 1 (x 2 )]=0 for (x l -x 2 ) 2 <0 (4.1) 

has causal support on the (closed) four-dimensional light-cones supp(D 2 ) = V + (0) UV~(0) (see 
Fig. 2), which are defined by 

V + (x) = {y I (y-x) 2 >0,y°> x }, V~ (x) = {y \ (y-x) 2 > 0, y° < x }. (4.2) 

(') (') 

One therefore introduces (primed) advanced and retarded distributions A 2 (z) and R 2 (z) according 
to 



R 2 = +D 2 A 2 

V+-{0] 



-D 2 _ , R 2 ,A 2 

v--{o} 

Ti(x 2 )Ti(x]) , A' 2 = -7i(xi)7i(x2). 



elsewhere, 



(4.3) 
(4.4) 



The non-trivial (!) splitting of D 2 into the retarded and advanced distributions 7? 2 and A 2 (see 
Fig. 2) indeed corresponds to time-ordering: One has T [T\ (xi ) T\ (x 2 )] " = " T 2 (xi , x 2 ) = R 2 — R' 2 = 
A 2 — A' 2 . One may substantiate this observation by explicitly checking that /? 2 — R' 2 is given for 
z° > by 7i(xi)7i(x 2 ) - ri(x 2 )7i(xi) + T l (x 2 )T l (xi), and for z° < by +T l (x 2 )T 1 (x l ). However, 
the decomposition is not unique in general in the critical point z = 0. 
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t z° 


v + ^ 






4- s» 




\ Z 1 ' 2 




supp(A2)j) 







Fig. 2: D 2 (z)=R 2 {z) -A 2 (z) =# 2 (z) -A' 2 (z) 
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Considering <p 3 -theory as an example with T\(x) = ^ : <p(x) 3 : and the scalar Feynman propagator 
defined via (0|r(<p(xi)<p(x 2 ))|0) = iAp (x\ —X2), standard Wick ordering leads to 

A 2 9A 2 
r 2 (xi,x 2 )" = " - ^ : <7>(*i) 3 <7>(*2) 3 : — : <p(xi) 2 <p(x 2 ) 2 : iA F (xi -x 2 ) 

18A 2 A 2 
— ^2" : 9(xi)(p(x2) : [iA F {xi -x 2 )] 2 - — [iA F {xi -x 2 )] 3 . (4.5) 

However, in the causal approach one constructs first 

A 2 9A 2 
D 2 (xi -x 2 ) = <p(x!) 3 :,: (p(x 2 ) 3 :] = ...- — : <p(xi) 2 <p(x 2 ) 2 : iA(xi -x 2 ) + ..., (4.6) 

where A(xi — x 2 ) is the Pauli-Jordan distribution, which can be decomposed into the positive- and 
negative-frequency Pauli-Jordan distributions A(z) = A + (z) + A~(z) given by 

A ± (z) = J d 4 ke(±k°)S(k 2 -m 2 )e- ikx . (4.7) 

4.2.1 Tree level 

In order to get the retarded C-number part of d' 2 ree {z) := A(z) of the particle -particle scattering 
diagram, we simply multiply by &(z°): r'{ ee (z) = &(z°)A(z). From &(k) = i ^ ) 8 { ^(k) follows 
(in a 'sloppy' style) in momentum space 



(4.8) 



2 vv ' " 2nJ 1 k°-p° + iO 2nJ 1-r+fl)' 
or, from the Lorentz covariance of A(z) (A(k)) one may derive a dispersion relation 

f"W-s/*T^i) fork€V+ ' (49) 

From A(k) = -2nisgn(k°)8(k 2 - m 2 ) follows 

J 2Vi?m(\-t + iO) k 2 -m 2 
The full expressions for r 2 ree (fc) and t 2 ee {k) for arbitrary k follow from analytic considerations. 

4.2.2 Loop level 

The self-energy (one-loop) part in T 2 is logarithmically divergent: 

t l 2 °° p {x l -x 2 ) ~ [iA F {x\ -x 2 )] 2 [ — 2 ,.J" P , 2 2—^. (4.11) 

z J [p 1 -m 1 + iQ\[{k- p) 1 -m z + iO\ 

In the causal approach, one calculates first D 2 (x\ — x 2 ) = [Ti(xi),7i(x 2 )], leading to 

d l 2 oop (x l -x 2 ) ~ [A-(xj -x 2 )] 2 - [A-(x 2 -xO] 2 sgn(£°)0(£ 2 -4m 2 ). (4.12) 
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Naive splitting leads to a divergent dispersion integral, as in the Feynman integral case: 

&(z°)d l r(z) ^ f dtf^- Q (4.13) 

However, it can be shown that the retarded part of d 2 ° op can be obtained from a finite, subtracted 
dispersion integral (in the massive case m / 0) 

*"<* ) = £/ (,- M )^(i-«+fl) + "'" '■ (iev+) wi,hffl=a <414) 

— CO 

(0 derives from the scaling properties (— ► power counting degree of divergence) of the causal dis- 
tribution d l 2° p . 

4.3 Higher orders 

With causality as the fundamental input, it is possible to construct causal distributions 
A n , R n and D n {x\ ,x n ) = R n (xi, . . . ,x n ) — A n {x\ ,x n ) at every order of perturbation theory, 
which have causal support 

supp/?„ (*i,...,*„) cr + (x„), suppA„(xi,...,x n ) C F~(x n ), (4.15) 

where the generalized closed forward (backward) light-cones are defined by 

r ± (x n ) = {(x u ...,x n )\x j eV ± (x n )Vj=\,...,n-\}. (4.16) 

Obviously, one has suppD„(xi, . . . ,x n ) C T + {x n ) ur(i„). We give here a short recipe how the 
higher order distributions are constructed. The crucial step in the inductive construction of the T n is 
the splitting of the D n , which can be performed in a mathematically well-defined ('finite') manner. 
However, the result is in general not unique, a fact which is directly related to the possibility of 
finite renormalizations of Green's functions in the standard approaches to renormalization theory. 
Recipe: T m known for 1 < m < n - 1 

construct primed distributions 

A n (x\, . . . ,x n ) = / T m (X)T n _ n , (Y,x n ) and R n (x\, . . . ,x n ) = } T n - n , (Y,x n )T m (X) 
p 2 p 2 

withP 2 : {jci,...,jc„_i} =XUY,X ^<d,ni = \X\ > 1 

allow X = 0, consider distributions with causal (light-cone) support 

An{%\ j ■ • ■ — A n [x\ , . . . + T n {x\ , . . . , R n (x\ , . . . ,x tl ) — R n [x\ , . . . ,x tl j + T n (x\ , . . . 
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T„ unknown, but difference distribution 
D n = R n — A n = R n — A n 
can be shown to be causal: distribution splitting of D n generates R n 

n> 

T„ = R n — R n ! 

As demonstrated above at second order for the particle -particle scattering tree diagram and the 
self-energy loop diagram, the C-number parts of the D n , R n : d tree ' loop, "(x\ — x n , . . . ,x n -\ — x n ), 

/r eeJoo P ,...( Xi _ %n ^ _ _ ?JCni _ %n) goover into ftree,looPr-( pu . . . jP ^) an d d tree ' loo P'-(p U . . . ,p„_i) 

via Fourier transformation. If at least one field in a field theory is massive, it can be shown in 
general that at all orders of perturbation theory, for p = (p\ ,p2, ■ ■ .) £ T + a subtracted dispersion 
relation applies as a distribution splitting operator 

*M = TnS (.-nr^L + nf +Z («^ulti-mdex), (4.17) 

— oo 1^1 ^ 

where ft) is a rigorously defined power-counting degree of divergence. The terms Y, c aP a corre- 
spond to possible finite renormalization terms (which have been discussed in connection with the 
renormalization group of QED in [3]). The polynomial terms in Fourier space are due to the fact 
that the splitting of D n into R n and A n is not uniquely defined at the "tip" x\ = X2 = . . . = x n of 
the generalized forward/backward light-cones r . In configuration space, they correspond to lo- 
cal terms ~ Y,CaD a 8(xi — x n , . . . ,x n -\ —x n ). Perturbation theory alone does not specify them, 
and they have to be restricted, e.g., by symmetry considerations and may lead to subsequent finite 
renormalizations of the T n , as already mentioned. 

5. Gauge theories 

We shortly mention the case of purely gluonic QCD as a standard example for a gauge theory. 
Within the causal approach, it is most natural to start from a first order gluon field coupling (matter 
fields neglected) 

Ti (x) = i\f ahc : A^(x)A b v (x)F c ^(x) :, F^x) = d»A v c (x) - d v A?(x). (5.1) 

First order gauge invariance requires additional fields (ghosts) 

Ti(x) = igfabc : A^(x)A h v (x)d v A^x) : -igf abc : A^{x)u b {x)d^u c {x) : . (5.2) 

At second order, tree diagrams containing C-number distributions ~ d^d v A(xi —X2) -^ — k^k v k{k) 
appear; their splitting is fixed up to a local term ~ g^ v 8^{x\ —X2) according to distribution theory. 

k^k v 

k^k v A{k) s P^ 4r^+Cg^ v . (5.3) 
k l + z'O 

Second order gauge invariance determines the constant C such that the usual four-gluon coupling 
term is automatically generated [12]. 
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6. Concluding remarks 

We conclude by giving a short and incomplete review of recent activities concerning the causal 
Epstein-Glaser approach. The causal approach was 'rediscovered' by Michael Diitsch and Giinter 
Scharf (U. Zurich) in 1985. Their work resulted in a textbook [3], where a complete discussion 
of QED and introduction to the causal approach can be found. From 1989-1993, the Zurich group 
led by Giinter Scharf focused on, e.g., interacting fields in the causal approach [4], axial anomalies 
[5] and gave a full discussion of the renormalizability of scalar QED [6]. During the period from 
1993-1999, a complete discussion of perturbative QCD was worked out [7], and gauge theories 
like the full standard model (including the phenomenon of spontaneous symmetry breaking) were 
studied [8, 9], theories in dimensions other than four were considered [10] and analytic calculations 
of multi-loop diagrams were carried out [11]. Quantum gravity and supersymmetric theories were 
also considered, and other groups (Klaus Fredenhagen et al.) generalized the causal approach to 
field theories on curved space-times. Finally, recent work of Ernst Werner and Pierre Grange in 
connection with light cone quantum field theory should be mentioned, which will also be part of 
the LC2008 proceedings (see also [13]). 
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